AD-A126  309 


UNCLASSIFIED 


SPACEFILLING  CURVES  AND  ROUTING  PROBLEMS  IN  THE  PLANE 
(U)  GEORGIA  INST  OF  TECH  ATLANTA  PRODUCTION  AND 
DISTRIBUTION  RESEARCH  CENTER  L  K  PLATZMAN  ET  AL. 

11  FEB  83  PDRC'83-02  N000 14-80-K-0709  F/G  12/1 


1/  \ 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NA1IONAI  BuRlAU  01  SIANOAROS  |96.<  A 


f 


1 


m  FILE  C0PY 


jID  A  X  2  6  3  0  9 


•<*  *i 


JO1/  t '  ■  '■ 


SPACEFILLING  CURVES  AND 
ROUTING  PROBLEMS  IN  THE  PLANE 


LOREN  K.  PLAT2MAN 
JOHN  J.  BARTHOLD1 1  III 


PDRC  93-02 


M 


:  ’  •v'- 
;  *  v\,  • 

•  -  ,*  j. 

V:<  •  -d^-f  v 


w  /  •  •  .  . 

\  V-*5*  .«• 

r  ■  >  X  •  T-  '• 

'  .  :>}'  V  :  :  ■ -  ..  ... 

DTIC 

*  tr»  cr  rt 


APE  4  1963 


PRODUCTION  and  ^  A 
DISTRIBUTION  RESEARCH 


Copt  available  to  DTIC  doe*  not 
fully  legible  reproduction 

(Thi*  document  has  been  approved 
ior  public  release  and  sale;  its 
distribution  is  unlimited. 


83  08  22  012 


•*i£>3r 


PDRC  Report  Series  83-02 
February  11,  1983 


Contract  N00014-80-K-0709 


SPACEFILLING  CURVES  AND 
ROUTING  PROBLEMS  IN  THE  PLANE 

LOREN  K.  PLAT2MAN 
JOHN  J.  BARTHOLDI  III 

PDRC  83-02 


O’hC 

■  -4  ELECT"  F 


APR  4  1983 


f 


School  of  Industrial  and  Systems  Engineering,  Georgia  Institute 
of  Technology,  Atlanta  GA  30332.  Telephone  <40 4 >-894-2300 

SUBMITTED  TO  MATHEMATICS  OF  OPERATIONS  RESEARCH 

This  document  was  prepared  on  a  word  processing  system  devised 
by  Loren  Platzman.  It  runs  on  the  Radio  Shack  Model  III  with  a 
very  cheap  printer.  Please  let  us  know  what  you  think  of  it  I 


This  document  has  bean  approved 
lot  public  telease  and  sale;  its 
distribution  is  unlimited 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DTIC  CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


SPACEFILLING  CURVES  AND  ROUTING  PROBLEMS  IN  THE  PLANE 
Loren  K.  Platzman  and  John  J.  Bartholdi  III 
Georgia  Institute  of  Technology 

This  paper  introduce*  a  novel  heuristic  to  compute  a 
minimal-length  tour  of  N  given  points  in  the  plane*  they  are 
sequenced  as  they  appear  along  a  spacefilling  curve.  The 
algorithm  consists  essentially  of  sorting*  It  is  easily  coded; 
requires  only  0(N)  memory;  and  may  be  implemented  to  execute  in 
0(N  log  N)  operations  at  most;  or  0(N)  operations  on  the 
average. *  If  the  points  lie  in  a  square  of  area  A;  the  heuristic 
tour  will Vave  length  2  JNft  at  most.  If  the  points  are 
statistically  independent  under  a  smooth  distribution;  with  N 
large;  then  the  tour  will  be  approximately  25*  longer  than 
optimum;  and  a  simple  enhancement  reduces  this  to  15H.  Me  also 
give  performance  bounds  for  our  method  whan  the  given  points  li* 
in  a  general  subset  of  d-space;  with  an  arbitrary  distance 
measure. 

AMS  1980  subject  classification.  Primary  *  90C42.  Secondary  * 

90033. 

QR/MS  Index  i960  subject  classification.  Primary*  491 
Network  s/*g  rap  hs/tr  aval  ling  salesman,  secondary  *  632 
Pr og  r amm i ng ✓ i nteg  er /a 1 g  or 1 thms^heur i Stic. 

Kau  words*  Planar  travelling  salesman  problem;  heuristic 
algorithm;  spacefilling  curve. 

The  authors  have  been  supported  by  the  New  Faculty  Research 
Development  Program  of  the  Georgia  Institute  of  Technology;  and 
bu  the  Office  of  Naval  Research  under  Contract  No. 

N000 14-80 -k -0709.  Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  U.S.  Government. 


1.  Introduction.  The  travelling  salesman  problem  <TSP) 
is  to  construct  a  circuit  of  minimum  total  length  that  visit* 
each  of  N  given  points.  Even  in  the  plana,  this  problem  is 
HP-complete  C123.  Karp  £103  has  given  an  CKN  log  N> 
heuristic  to  compute  a  tour  whose  length  is  within  i  of 
optimal,  but  it  is  difficult  to  code,  and  its  effort  has  a 
constant  factor  that  Increases  rapidly  as  «  decreases.  Bentley 

g/o 

and  Saxe  C63  devised  an  CKN  log  N)  implementation  of 
the  nearest  neighbor  heuristic,  but  it  requires  a  special  data 
structure.  This  paper  describes  a  faster,  simpler  heuristic 
that  performs  comparably.  He  gave  a  brief  account  of  the  ideas 
underlying  our  method  in  C33.  Here  we  present  the  algorithm 
in  detail  and  analyze  its  performance. 

Let  U  be  a  set  <e.g..  the  unit  square?  within  which  tours 
are  to  be  constructed,  and  define 

C  *  <6  I  0  <  8  <  1>.  a.l) 

When  U  has  dimension  2  or  more,  a  continuous  mapping  *  from  C 
onto  U  is  known  as  a  spacefilling  curve.  Such  curves  were 
first  devised  by  Peano  and  Hilbert  in  the  1890's  to  resolve 
topological  existence  questions  <cf.  Hobson  C9.  pp  451-4583). 
More  recently,  they  have  attracted  attention  as  entertaining 
examples  of  recursively  defined  computational  procedures  Ccf. 
fllaph  Naught C23.  flbelson  and  diSessaCl.  pp  94-1023).  We  turn 
them  to  decidedly  practical  use. 


Let  v<0>  *  1  > ,  Then  'KB)  trace*  out  a  "tour"  of  all 

the  point*  in  U  at  9  varies  from  0  to  1.  CFor  this  reason/  we 
find  it  convenient  to  view  each  6eC  at  a  point  on  the  unit 
circle/  9  clockwise  revolutions  from  a  fixed  reference  point. > 
Given  N  points  in  u  to  be  visited/  our  strategy  is  to  sequence 
then  as  they  appear  along  the  spacefilling  curve.  Thus/  we 
propose  to  construct  tours  in  the  following  Manner < 

SPACEFILLING  HEURISTIC 

1>  For  each  point  $  to  be  visited/  compute  a  0  such  that  $»>K8>. 
2>  Sort  the  points  by  their  corresponding  0's. 

In  Section  2/  we  analyze  the  performance  of  this  heuristic 
when  the  region  U/  the  distance  measure  by  which  tours  are 
evaluated/  and  the  spacefilling  curve  V  are  all  arbitrary. 

Then/  in  Section  3/  we  define  a  particular  spacefilling  curve 
over  the  unit  «nuare/  and  show  how  its  inverse  is  computed.  This 
provides  a  specific  algorithm  to  solve  TSP's  in  the  plan®;  its 
performance  is  discussed  in  Section  4.  Section  5  shows  how  the 
ideas  of  Section  4  can  be  extended  to  more  general  TSR  problems. 
Sections  6-7  contain  lengthy  proofs  deferred  from  earlier 
sections.  Concluding  remarks  are  given  in  Section  3. 
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2.  General  performance  analysis.  Suppose  that  U  lit* 
in  d-space  and  1st  11*11  b«  a  norm  that  da  tor  Mi  no*  the  distance 
between  points  in  U.  Our  analysis  is  based  on  the  following 
three  technical  conditions,  which,  as  we  shall  see,  hold  for 
most  problems  of  practical  interest1 

CPI >  The  inverse  of  v  is  easily  evaluated.  Specifically, 
if  t«U  is  a  d- vector  whose  components  each  have  k-bit 
representations,  then  there  is  a  8  satisfying  **+ <8 >  and 
having  an  0(dk>-bit  representation  which  may  be  computed  in 
0<dk>  operations.  Although  there  may  be  many  such  8,  we  must 
compute  only  one. 

<P2)  There  is  a  finite  constant  ii  such  that 

llv<0)-t<8 '  >U  <  f<l8-6'l>,  8,6'  eC,  <2.1> 

where 

f<A>  -  <1  CminCa,  i-a>J1/'d,  a«C,  <2.2> 

<P3>  *  is  Lebesgue  measure  preserving.  That  is,  if  I  is 
an  interval  in  C,  the  tet  <  +<8>  I  0«I  >  has  d-volume 
(area  if  d»2>  equal  to  the  length  of  I. 


4  ~ 


Remarks.  <P1>  ensures  the  efficiency  of  our  algorithm/ 

<P2)  enables  us  to  study  its  performance;  and  <P3>  normalizes 
the  problem  and  eliminates  patholog ical  cases.  Intuitively / 

<P1>  »$■  t  is  surjective  and  <P2>  +  is 

continuous/  so  together  they  imply  that  ♦  it  i  spacefilling 
curve.  Also/  <P2)  v<0 >  *  <Kl>.  <P3>  requires  that  U 

and  C  be  scaled  so  that  U  has  unit  d~volume  and  * 

"homogeneously  fills"  U. 

Computational  Effort. 

The  spacefilling  heuristic  requires  0<dkH>  operations  to 
compute  the  <?'s  and  0<N  log  N>  operations  to  perform  the  sort/ 
that  is.  OCN  log  N)  operations  in  all.  (To  be  more  precise/  dkN 
is  the  number  of  bits  required  to  specify  the  probleM/  and  the  N 
numbers  to  be  sorted  each  have  0<dk)-bit  representations/  so  the 
problem  is  solved  in  OCdkN  log  N>  bit  operations.) 

Alternatively,  a  randomized  sorting  algorithm  requiring  0(N> 
expected  comparisons  might  be  used  <e.g.  B INSORT/  see 
KnuthCil]  and  WeideC153).  A  heuristic  tour  is  then  obtained 
in  O(dkN)  expected  operations. 

This  procedure  also  requires  surprisingly  little  memory/ 
only  0<dN>/  linear  in  the  size  of  the  problem. 
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Worst-case  ftnalusis. 

Given  N  sorted  points  +<8^*  •  *>K8^>,  the  tour  length  is 
bounded  above  by  where  Ai-8i+i~®i'  i“i,***,N-l 

and  Since  this  expression  is  concave  and  symmetri 

_,u 

in  and  since  y^j  Ai  *  1,  it  achieves  a  maximum  of 

NfCi/N).  Thus,  we  have  proved 

Theorem  2.1.  The  spacefilling  heuristic  produces  a  tour 
whose  length,  L,  satisfies 

L  <  Nfa^hO  -  0  M<d"1>/d. 

•m 

We  also  consider  how  long  the  heuristic  tour  nay  be  in 
relation  to  the  optimal  tour. 

Theorem  £.2.  If,  for  a  given  set  of  N  points,  the  heuristic 

£ 

tour  has  length  L  and  the  optimal  tour  has  length  L  ,  then 
L-'L*  <  0<log  N>. 

The  proof  of  Theorem  2.2  is  lengthy  and  will  be  given  in 
Section  8.  We  note  that  the  dominant  term  in  the  exact  bound 
for  LrL*  <Eq.  <6. 5))  increases  slowly  in  N  but  rapidly  in  d. 
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Probabilistic  Rnalysis. 

Suppose  that  the  points  *j*  ^  are  independent 

random  variables  uniformly  distributed  over  U.  Mow  the  heuristic 
tour  length  and  the  optimal  tour  length  may  be  viewed  as  random 
variables.  Since  <¥  is  measure  preserving*  each  8i  satisfying 
is  uniquely  determined  with  probability  one* 
and  is  uniformly  distributed  over  C.  When  N  is  larger  the 
6^'s  approximate  a  Poisson  process  on  C. 

Define  random  variables 


n* 

PN 


N-<d-l'j/d  ,  # 


N» i  >  2^  •  1  ■  ; 


where  L*  is  the  length  of  an  optimal  tour  joining  M  random 
points.  Beardwoed*  Hal ton  and  Hammer* l ay c 5 3  have  shown  that- 
when  the  Euclidean  norm  is  used  to  evaluate  tours*  there  is  a 
p*  such  that  lim^^pjjj)  ■  p*  >  almost  surely.  Steele  C 143  has 
extended  this  to  the  stronger  statement  of  convergence 


£n»1  ProbCifH  "  «  «>0. 


Me  now  establish  a  similar  result  for  the  heuristic  tour. 
Define  random  variables 


pH 


tj**  C  d—  1 )  rd  1  lj_<  n 

N  L  *  N* 1 j  £*...* 


<2. 3) 


where  L  is  the  heuristic  tour  length  for  M  random  points*  and  let 


pN  -  ECpN3. 


<2. 4> 
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The  sequence  need  not  converge  a*  N*«*.  By  Theorem  2.1. 
however.  it  is  bounded1 

0  <  <  A.  <2.3> 

These  bounds  may  be  tightened  slightly  by  noting  that  (i>  the 
optimal  tour  grows  as  n^-lVd  p* >  so 

p*  <  11m  inf^  <M-^> ■*  <2  - 

and  <ii>  the  increments  in  $  along  a  heuristic  tour  are  nearly 
exponentially  distributed  with  mean  1/N.  so.  by  concavity  of 
fCO. 


lim  supM  <pN> 

<  h  J  f<x>  Ne“N*  dx 

■  A  r<<d-lVd).  <2.7> 


Therefore.  for  large  N.  the  expected  heuristic  tour  length  will 
be  within  a  constant  factor  A  r<Cd-l)/d>  /  p*  of  the 
expected  optimal  tour  length. 
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The  following  theorem  will  permit  us  to  m«.k *  a  stronger 
statement  rag  ar  ding  the  ratio  of  Kauri  at  ic  to  optimal  tour 
lengths. 

Theorem  2.3.  Lat  ba  a  (deterministic)  sequence  such 
that  ■*  0.  Than 

*  *■*  *•'  W  «>0. 

Proof »  Ua  slightly  modify  Steele's  proof  C141  of  complete 
■* 

convergence  of  pjj.  Simply  radafina  the  random  variables 
d^  in  C143  to  reflect  the  fact  that  our  heuristic  tour  joins 
a  point  with  another  whose  inverse  image  under  f  is  closest  in 
C  under  metric  f<0; 

dN  ■  nin<f<l0i-81l>  I  im2t . . . >N>. 

Since  the  8^'s  are  independent  and  uniformly  distributed  on 
C>  Eq.  <2. 2)  of  C143  may  be  replaced  by 

ProbtdN  >  t3  <  Cl-2f"lCt>3N"1  <  Cl-2<t/ft>d3N“1/  0<t<fCl/2>. 

and  the  remainder  of  C143  establishes  the  desired  result.  □ 


~  9  ~ 


*  * 

Sine*  p^  -  and  pj^  -  p  both  converge 
comp  lately  to  zero  as  N  ■*  *>  the  difference  between 

nil 

pN/pN  random  ratio  of  heuristic  to  optimal  tour 
length*)  end  <the  deterministic  ratio  of  expected 

heuristic  tour  length  to  expected  optimal  tour  length)  vanishes 
Almost  surely  as  N*».  Thus  the  spacefilling  heuristic  will 
reliably  construct  tours  whose  length  is  a  fixed  constant  factor 

nK 

l*r9#r  than  optimal.  This  factor  depends  on  the 

region  U>  the  curve  v.  the  distance  measure  IMi.  and  (in 

principle  although  not  in  our  experience)  the  problem  size  N. 

We  can  estimate  by  Monte-Carlo  simulation*  generate  many 

random  problems  and  solve  them  by  the  spacefilling  heuristic. 

£ 

Good  approximations  to  p  are  given  in  C53  for  problems  in 
subsets  of  d-space  under  the  Euclidean  norm.  For  d"2. 
p *  «  .765  C63. 
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Rn  addition*!  performance  measure  concern*  the  length  of 
tha  longest  link  along  a  tour.  In  general  tour*/  it  may  be  close 
to  tha  dianatar  of  U.  This  can  laa  to  difficult!**  in  certain 
application*.'  and  ha*  motivated  tha  formulation  of  the 
bottleneck  TSP  problem  C73«  to  construct  a  tour  who** 
largest  link  ha*  minimal  length.  Our  hauri*tic  produce*  tour* 
who**  longest  link  i*  not  much  greater  than  tha  average  link 
length*  a*  we  now  «how. 

Theorem  2.4.  The  longest  link  along  a  spacefilling 
heuristic  tour  has  expected  length  0<log<‘.M>''N>1//d. 

Specifically;  if  t>  1*  the  length  of  the  longest  link; 
lim  *upN4w  E((ln(N)/N)'l/d  6)  <  ft. 

Proof  <  Let  a  be  the  largest  increment  in  8.  The 
increment*  in  8  may  be  considered  a*  Independent  exponential 
random  variable*  of  mean  1/N;  *o  the  probability  that  none  of 
the  N  increment*  exceed*  t  is  <l-*xpONt>>N.  Let  t*ln<N</sVN. 

The  probability  that  all  N  increment*  are  bounded  by  t  become* 
<l-*/tON.  fl*  H  increa*e*  thi*  become*  expo*).  It  follow* 
that  a  i*  distributed  a*  ln<N)XN  ~  in<*)xN;  where  •  i*  an 
exponential  random  variable  of  mean  1.  So  a*  N+*; 

E<A*'nn<N>'/N3>  ♦  1.  Since  f  is  concav*; 

E<&>  <  E<f<A>>  <  f<E<4>>.  Q 


A 
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3.  R  Spacefilling  Curve  in  the  Unit  Square.  We  now 
construct  a  curve  that  fills  the  unit  square* 

S  -  <  <x*y>  I  0<x<l,  0<y<  1  ).  0.1) 

It  resembles  Sierp inski's  curve  c 133  and  enables  our  heuristic 
to  perform  especially  well.  It  is  recursively  defined  by 
breaking  8  into  four  identical  subsquares*  and  filling  each  with 
a  spacefilling  curve  rotated  so  that  the  four  subcurves  link  to 
form  a  circuit  (Figure  1). 

We  may  arbitrarily  specify  the  "starting  point"  of  the 
curve  to  be  <0*0)*  and  the  "direction"  of  the  curve  to  be 
clockwise.  Then*  as  0  increases  from  0  to  l*  the  four 
subsquares  are  visited  in  the  sequence 


so 

■ 

< 

(x*y)cS 

1  x<.3* 

y<.5 

> 

S1 

m 

( 

<x*y >«S 

l  x<.  5* 

y>.3 

> 

S2 

■ 

< 

(x.yJtS 

1  x>.  5* 

y  >.  5 

> 

S3 

■ 

< 

<x*y)«S 

1  x>.  5* 

y<.5 

>. 

<3. 2) 
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These  subsquares  induce  a  partition  of  C  into 


C0  *  <  0  i  0<8£i/e  or  7/8<8<i  > 

Cj  ■  <  8  i  l/’8<8<3/,8  > 

Cg  -  <  8  I  3/8<8<3/8  > 

Cg  *  <  8  I  5/8<8<7/8  ). 


0.3) 


so  that  4  naps  CA  onto  St*  i  ■  0*i*2>3.  For  ease  of 
definition  we  need  an  operation  representing  the  rotation  of 
each  subsquare.  Let  R^  be  an  affine  transformation  from  8  to 
8^  determined  by 


<0, 

.0) 

♦ 

<.3* 

.  5>* 

<0* 

1> 

4  <.5i 

.0)* 

a 

*1> 

♦  <0, 

-0>* 

<1 

*0)  4 

<0* , 

.3) 

<0. 

.o> 

4 

<.3* 

.s>. 

<0* 

1> 

4  <0  *  , 

,3>* 

<1 

*  n 

4  <0; 

.  1>* 

a 

*0>  4 

< .  5j 

.  1? 

<0. 

,0) 

4 

<.3* 

.  3)  * 

<0* 

i) 

4  <.5, 

1>* 

<1 

*  i) 

♦  <1. 

.  1>* 

<i 

i0>  4 

<1* , 

,5> 

<0, 

.0) 

4 

<.3* 

•  5>* 

(0* 

i> 

4  a,. 

,3>* 

<1 

*  1) 

4  <1, 

,Q>* 

a 

*0>  4 

<.5* 

,0) 

(3.4> 


These  are  interpreted  as  follows1  position  <x*y>  relative  to 
subsquare  i  corresponds  to  position  R^sy)  relative  to  the 
larger  square  8.  For  exawple*  the  point  <0*0  >  marks  the  start 
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of  tha  spacaf illing  curves  to  ^<0,0)  is  tht  starting  point  of 
tha  aubcurva  ovar  S . ,  that  is,  tha  eantar  of  S»  <.3 ,.3). 
Finally,  dafina  afftna  transformations  from  aach  ci  onto  C' 


qQ  •  7/8  0,  1  * 
qj  '  1/8  ♦  0,  1/4 
q2  !  3/8  ♦  0,  1/2 
q3  1  5/8  ♦  0,  3/4 


.5;  0  *  .5,  1/8  ♦ 

♦  .3,  3/8  ♦  1 

♦  .5,  1/8  +  1 

♦  .5,  7/8  +  1. 


1 

0.5) 


Tha  purposa  of  qA  is  to  convart  positions  along  tha  largar 
curva  <ovar  8)  into  poaltlona  along  tha  aubcurva  ovar  8 . .  For 
axanpla,  Sj  ia  covarad  by  •¥  aa  0  rangaa  ovar  C^,  and  0».3 
ia  ancountarad  qjC.3)*.?  of  tha  way  through  tha  aubcurva  ovar 
Sj.  Using  this  notation,  tha  racuraiva  atructura  of  <v,  shown 
in  Figura  1,  is  concisaly  axprasaed  as 


*<e>  -  RiCt<qiC03)3,  eiCj, 


i*0, 1,2,3.  0.6) 
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Figura  1  auggeat*  that  a  hauriatic  tour  will  conaiat 

aaaantially  of  four  aubtoura,  aaeh  involving  a  fourth  a*  Many 

point*  in  a  fourth  a*  larga  an  araa.  Conaaquantly,  ua 

anti cl pats  that  tha  axpactad  hauriatic  tour  langth  for  4N  random 

point*  in  8  will  ba  twica  a*  long  a*  that  for  N  point*.  That  ia, 

**4N  "  ^N>  ***  *"d  <2.4).  Vat  taquanca*  of  tha 

form  <p.  u>*  might  approach  a  limit  £(•>  that  dapand* 

N*4k  k»l 


log4<N>  mod  1 

on  Ni  or,  Mora  *p«cif ical ly ,  on  4  * 


Such  i* 


tha  caaa,  a*  wa  aha 11  «aa. 


It  follow*  from  0.6)  that  •¥>  if  it  axiata,  aatiafia* 


♦  1  0  <0,0>,  1/4  (0,1), 


1/2  ♦  (1,1).  3/4  ■>  (1,0), 


(3.7) 


Eq*.  0.6)  and  0.7)  provida  tha  ba*i*  for  tha  following 
algorithm  to  conputa  tha  invar*a  of  «♦" 


FUNCTION  THETR<X,Y) 

2.  Chooaa  I  (not  nacaaaarily  uniqua)  so  that  <X,Y)iSj 
2.  If  X«<0,1>  and  Yi<Q,i> 
than  RETURN < I /4> 

Othar w i * a  RETURN Cq j 1 C THETfl <  R  j 1 C  X ,  Y  3  )  3  ) . 
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The  following  theorem  summarizes  proper tie#  of  *  that  may 
reasonably  be  inferred  from  our  discussion  thu#  far.  It*  formal 
proof  1*  lengthy,  however,  end  will  be  given  in  Section  e. 

Theorem  3. 1 .  The  function  v  defined  by  <3. 6)  exist* 
uniquely  And  i*  spacefilling  on  S.  It  satisfies  <Pi>  (its 
inverse  being  computed  by  the  function  THETfi),  <P2>  with  IMI 
»  Euclidean  distance  in  (2.1)  and  A  »  2  in  (2.2),  and  <P3>. 
Moreover,  there  is  a  continuous  function  p  on  Cl, 4)  such  that 

,  log^tN)  mod  1* 

iiP,N-r«  <f*N  ~  p^4  J>  -  0, 

when  corresponds  to  tours  in  6  based  on  the  curve  defined 
by  0.6). 
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4.  Performance  assessment  for  problems  In  the  square. 

To  obtain  a  Kauri at lc  tour  in  a  square  of  araa  R.  we  raacala  tha 
pointa  ao  that  they  11a  in  tha  unit  square.  compute 
corresponding  8  values  by  tha  procedure  THETR  of  aaction  3*  and 
aort.  <fl  abort  BASIC  coda  for  thia  algorithm  was  giv/an  in 
C31«>  Tha  tour  obtained  ia /  of  course/  4TT  timaa 
long ar  than  a  tour  In  tha  unit  aquara.  Wa  now  compara  ita 
performanca  with  that  of  tha  naaraat  neighbor  heuriatlc  <NNH> 
and  tha  minimal  apanning  traa  hauriatic  <MSTH> .  aa  analyzed  by 
Bant lay  and  Saxe  C63. 

Computational  effort. 

Tha  worat-caaa  effort  to  generate  apacafilling  heuriatlc 
toura  ia  0<N  log  N)/  considerably  laaa  than  those  for  the  NNH 
and  MSTH-  CKN3^2  log  N)  and  0<N2  log  N>,  respectively. 

The  avaraga  affort  for  tha  NNH  using  a  special  data  structure, 
ia  0<N3/3),  but  tha  avaraga  affort  for  tha  apacafilling 
hauriatic  (using  BIN80RT)  la  0<N>.  Tha  computer  coda  for  tha 
apacafilling  hauriatic  ia  also  simpler/  and  requires  laaa 
memory. 

Morst-cass  performance. 

Both  tha  apacafilling  hauriatic  and  tha  NNH  produce  toura 
whose  length  is  at  most  0<-fT7FT>.  8o/  as  noted  in 
t63z  tha  ratio  of  hauriatic  to  optimal  tour  lengths  will  ba 
large  only  when  the  optimal  tour  ia  unusually  small.  Tha  worst 
case  of  this  ratio  la  O(log  N>  for  tha  apacafilling  heuriatlc 
with  general  U  and  ♦.  For  problems  in  tha  square/  however/  wa 
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suspact  that  It  is  bounded  bg  a  constant.  Th«  worst  ins tanca  wa 
hava  bun  abla  to  discover  Is  4.707;  for  a  32-point  problam.  fl 

L 

slmplar  axampla  Is  for mad  by  2  points  whosa  9 ' s  art  tha 
intagar  multiplas  of  2  j  tha  ratio  of  hauristic  to  optimal 
tour  langths  for  thasa  points  is  4.  Tha  worst-casa  ratio  for  tha 
NNH  is  known  to  lia  batwaan  0<log  N  /  log  log  N)  and  0<log  N)> 
for  tha  MSTH;  it  is  2. 

flvaraoe  performance. 

Tha  expected  length  of  a  hauristic  tour  of  N  points  is  N 
tlrias  tha  axpactad  langth  of  a  sing  la  link  along  that  tour.  Me 
ganaratad  10; 000  random  links  for  aach  of  various  valuas  of  N. 
From  thasa  experiments;  wa  concluded  that  fit ■ )  is  naarly 
constant  and  (probably)  lias  in  tha  ranga  .956  +  .05.  so 
that  tha  hauristic  tour  is  ".956  4 MW  .  Thus  our 
hauristic  producas  tours  "25*1  ovar  optimum,  fl  simp  la 
enhancement;  to  ba  givan  in  Sactlon  5.  raducas  this  to  15k  ovar 
optimum.  Tha  tour  obtainad  oy  tha  M8TH  is  ".95  /WT ; 

(  "25k  ovar  optimum);  and  that  obtainad  by  tha  NNH  is 
".92  'TRFT  ;  <  "20k  ovar  optimum). 
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Our  experiments  showed  that.  for  tours  produced  by  the 
spacefilling  heuristic,  the  longest  link  is  typically  between 
1.1  4  ?R/’HrTiS"R  and  1.3  4  (H/N)  in  N  .  The  constant  preceding 
the  radical  In  these  expressions  does  not  converge  as  N  ♦  ». 
Following  arguments  given  in  C73.  it  can  be  shown  that  the 
Minimum  possible  value  of  the  longest  link  has  expectation 
bounded  below  by  4  (ff/iWJU  TrT “TT  .  Thus 

our  method  provides  a  reasonable  solution  (within  an  average 
factor  of  2.3)  to  the  bottleneck  TSP.  Stated  differently,  the 
expected  ratio  of  longest  link  to  average  link  along  a  space¬ 
filling  heuristic  tour  is  0<4th  N  >.  The  length 
of  the  longest  link  along  a  NNH  tour  is  comparable  to  the  side 
of  the  containing  square,  so  the  expected  ratio  of  longest  link 
to  average  link  along  a  NNH  tour  is  0(4TT>.  (This 
ratio  is  significant-  in  applications  where  a  single  tour  is 
partitioned  into  subtours  for  each  of  K  vehicles.  I4e  would  like 
subtours  of  equal  length  to  contain  equal  numbers  of  customers. > 


-  20 


5.  Generalizations.  The  spacef ill Itiq  heuristic*  as 
descr ibed  In  Section  2,  it  not  restricted  to  the  method*  of 
Section*  3  «md  4.  Some  alternative*  ere  listed  below. 

Different  probability  distribution*. 

If  the  point*  are  not  uniformly  di*tributed  over  the 
square*  but  are  independent  with  identical  probability  densities 
n<x*y>  (and  N  i*  large)*  then  the  JTT term  in  the 
exp re** ion*  for  both  the  expected  heuristic  tour  length  and  the 

1  yO 

optimal  tour  length  i*  replaced  by  tj  <x*y)dy  dx 

(a  justification  l*  given  in  C53>.  Thu*,  the  ratio  of 
expected  heuristic  to  optimal  tour  length*  remains  unchanged. 

Other  metric*. 

Me  may  take  a*  our  measure  of  distance  (in  U)  the  1K  norm 

P 

H-i'flp  *  (  Zi-i  <*r*i>p  }l  P  P>* 

Thi*  include*  a*  special  case*  the  "sum  of  coordinates"  or 

"rectilinear"  metric  <p*l>  and  "largest  cooordinate"  metric 

<p«e) *  as  well  as  the  Euclidean  metric  (p*2).  Since 

HxIL  /  llxlU  lies  between  1  and  d*  1/2:> ,  * 

P  2 

curve  satisfying  <P2)  for  the  Euclidean  norm  also  satisfies  <P2) 
for  any  lp  norn. 
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Performance  enhancament. 

Our  Method  can  be  extended  * lightly  so  that  it  generate* 
better  tour*  in  the  square.  still  in  0<N  log  N>  operation*.  The 
additional  *tep.  briefly.  i*  a*  follow*,  After  computing  the 
•pacefilling  heuristic  tour.  attempt  to  move  each  point  to  a  new 
location  within  the  tour  by  computing  three  alternative  value* 
of  8.  More  specifically.  for  each  point  a  find 

the  largest  triangle  of  the  form  given  by  Lemma  7\5..  such  that 
lie*  in  the  triangle,  but  the  other  given  point*  do  not 
<l.e..  the  0  value*  for  the  remaining  point*  lie  outside  the 
interval  Ci»2  n.  <i‘*-l)*2-M3).  project  the  given 
point  onto  each  of  the  three  aide*  of  thi*  triangle  *o  that  the 
projection  lie*  Ju*t  outside  the  triangle,  determine  the 
po*ition  in  the  tour  of  each  projected  point,  and  move  the  given 
point  to  an  alternative  position  if  the  tour  would  be  *hortened 
by  doing  *o  (Figure  2?.  Me  have  estimated  the  length  of  the  tour 
obtained  by  thi*  improved  heuristic  to  be  ".88  .  that  is. 

"155*  over  optimum,  and  the  effort  to  obtain  an  improved  tour  to 
be  4-8  time*  greater  than  that  to  obtain  a  simple  spacefilling 
heuristic  tour. 


i 


When  a  point  is  shifted  onto  the  boundary  of  this  triangle,  its  position  in  the  tour  may  be  altered 
as  shown  in  the  second  square.  Returning  the  point  to  its  original  location  while  preserving  the 
new  sequence  yields  a  new  tour  of  the  original  points.  Retaining  only  those  sequences  that 
reduce  the  tour  length  produces  the  tour  shown  in  the  third  square.  Some  obvious  inefficiencies 
have  been  eliminated.  The  length  was  reduced  by  8/. 
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Other  spacefilling  curvt>  In  the  plane. 

Good  heuristic  tours  can  be  generated  from  spaeef i 1  ling 
curves  other  than  the  one  given  in  Section  2.  Hilbert's  eurve, 
for  example,  does  not  begin  end  end  at  the  sene  point;  this 
night  be  a  desirAble  feAture  in  sone  Applicstions  (such  as 
routing  a  bus  to  start  At  the  Airport  And  finish  downtown ,  then 
reposting  the  trip  in  reverse  order?.  If  Hilbert's  curve  is 
Arranged  so  that  it  starts  and  ends  at  the  same  point;  then  our 
experiments  show  that  it  will  yield  tours  assymptotical ly 
comparable  to  those  produced  by  ours.  For  smaller  N,  and  points 
distributed  over  a  square,  the  tours  constructed  by  our  curve 
are  slightly  better.  For  smaller  N  and  points  in  rectangular 
regions  (width  >  2* length?,  Hilbert's  curve  may  perform  better 
thAn  ours. 

ft  spacefilling  curve  in  d-space. 

0  curve  similar  to  that  of  Section  2  can  be  defined  by 
partitioning  the  unit  d-cube  into  2  similar  subcubas.  It 
visits  the  vertices  of  the  cube  (and  hence,  the  subcubes?  in  a 
sequence  determined  by  the  Gray  code  CS1 :  (0...0),  (0...0  1?, 
(0...0  1  1?,  <0. . . 0  1  0?,  (0...0  1  1  0?,  ...  .  If  we  take 
(0...0?  to  be  the  0-th  (and  2d-th?  vertex  visited,  then 
(1...1?  will  be  the  M^-th  vertex  visited,  where 


<2/3?*(2d-l?,  d  even 
(2/3?* <2d-.S?,  d  odd 
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Thus.-  for  d>2.-  <1...1>  will  not  be  visited  halfway  along  the 
curve.  Rather.  ■  <!...!>  for 


«d 


<2/3>*  a-2~al>.  d  even 
2/3.  d  odd 


Briefly.  ♦  visits  the  subcube  containing  the  n-th  vertex  as  8 
ranges  over  the  interval  C<n-ad>-2  .  <n+l-ctd>*2"’°3. 

n»l.- . . .  .2^-1;  and  Y<$>*the  n-th  vertex  when 
n  even,  or  8"<n-«d) *2“^.  n  odd. 

The  inverse  triage  under  ♦is  computed  by 


FUNCTION  THETfl<X1 ..... Xd> 

i>  For  i-1  to  d.  set  Vi  ♦  INT^fcX^.  Choose 
Nc<l. . . .  .2°*)  so  that  <Vj...VcJ)  is  the  N-th 
vertex  visited  by  the  curve. 

2>  For  i*l  to  d.  set  Y^  «■  1  -  2#min<Xi . 1-X^  > . 
CoMpute  8  «■  THETRCY^ ..... Y^> . 

If  N  is  odd/  then  set  8*1-8. 

RETURN <<N+e-ad>#2“d  mod  1>. 

Thus/  it  is  clear  that  <Rl>  is  satisfied. 
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To  show  that  <P2)  holds  <with  II *11  ■  Euclidean  distance), 
let  8j  and  8g  be  any  two  points  in  C.  and  set 
4  *  » 8 ^  — ® 2 1 «  D  »  ii'K81)-H'f.8g)ii.  Suppose  that 
the  cube  Is  broken  into  2kd  subcubes  of  side  2“k.  Me 
call  any  two  subcubes  neighbors  if  they  have  at  least  one 
vertex  in  common.  (By  this  definition,  a  subcube  is  its  own 
neighbor.)  Our  curve  visits  each  subcube  twice;  each  visit 
corresponds  to  an  interval  in  C  of  length  ctd  2~kd  or 
<l-ad>  2  and  consecutive  intervals  in  C  are  always 
mapped  into  neighboring  subcubes.  Thus,  if  a  <  min<ad.  1-cc^)  2‘"kd, 
then  y<8j)  and  +<8o>  lie  in  neighboring  subcubes,  and  so 
D  <  2  d1/£  2~k.  <Note  that  this  argument  is  valid 
only  for  integer  k.)  Thus.  <P2>  holds  with  ft  *  A  diYZ  <l-acJ)“1/'cl. 
<We  note  that  the  0  obtained  in  this  manner  is  not  as  good  as 
the  ft  established  in  Theorem  3.1  by  an  argument  that  is  valid 
only  for  d*2. )  Consequently,  the  spacefilling  heuristic  will 
produce  a  tour  of  length  at  most  OCd1^2  N<'d“1'’/d). 

Beardwood.  Halton.  and  Hammersley  C53  have  shown  that  the 
average  optimal  tour  is  of  this  order  of  magnitude. 
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6.  Proof  of  Theorem  2.2.  The  proof  is  based  upon  the 
following  identity1  if  A  is  a  sat  of  nonnegative  nunbars  and 
H<t>  *  #<  x«A  I  \>t  >/  than  tha  sum  of  alaments  in  A 
equals  /q  H<t)  dt.  We  let  A  represent  the  set  of  N  link 
lengths  along  the  heuristic  tour  and  construct  a  bound  of  the 
form  H*<t>  >  Wt>/  so  that 

L  ■  JC 

Let  S  denote  the  sat  of  N  given  points/  and  let  ECO 
denote  the  set  of  points  in  U  that  lie  within  distance  e  of  a 
point  in  S.  We  first  show  that 

V[E<<>3  <  c '  L*  *d"”^  +  c  td  <6. 1  > 

where  VC- 3  denotes  d-voluma/  c  *  VC<  %  I  II^IKl  >3  and 
cJ  »  supj^j^j  VC <  v’V+a'l  I  E^I^O/  0<a<l  >3. 

<If  IMI  is  Euclidean  distance.*  then  c  is  the  volume  of  a  unit 
d-sphare  and  c'  is  the  volume  of  a  unit  d-cylindir.  Note/ 
however/  that  II -  n  can  be  any  norm  on  d-spacei>  To  simplify 
notation/  let  "node"  indicate  a  point  in  S/  and  let  "link’1 
indicate  a  segment  of  the  optimal  tour  of  £.  R  point  will  be 
considered  to  "lie  within  «  of  a  link"  if  it  lies  within 
distance  <  of  some  point  along  the  link.  Since  any  given  set 
of  1  nodes  <i*l.* . . . /N-l>  adjoins  at  least  1+1  links/ 


*  H< t )  dt  <  H+Ct>  dt. 
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<%,  I  *  lit*  within  «  of  i  or  nor*  nod**) 
c  l  $  1 la*  within  <  of  i+1  or  wor*  link*)>  i«l,...*N-l 

It  follow*  that 

VE<%  I  %  lia*  within  «  of  1  or  nor*  node*) ] 

+  ^i*«i  VC<$  I  4  11a*  within  i  of  i  or  nor*  nod**) 3 

<  VC <4  I  %  lias  within  i  of  1  or  nora  link*)] 

+  VC<^  f  t.  li«*  within  <  of  all  N  nod**>3  <6.1'> 


We  nay  obtain  <6. 1)  fron  <6. 1')  by  tha  thr*a  *ub*titution*  yivan 
balow 

VC<*  I  \  lie*  within  i  of  1  or  nor*  noda»>3 

*  ^isl  i  •  VC<*  I  ^  lia*  within  i  of  i  noda*  exactly >3 

*  VC (4  I  *  lia*  within  «  of  tha  i-th  noda) 3 
-  N  C  «d 


1 
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E*  1  VC !  *  1  its  within  i  of  i  or  nor*  link*) 3 

m  j;N  i  •  VC<*  1  *  li**  within  «  of  1  link*  *xactly>3 
m  VC«  I  ^  li**  within  t  of  th*  i-th  l ink) 3 

-  c*  L*  «<d"1>  +  N  c  «d 

VC<*  t  %  li**  within  «  of  *11  N  nod** >3  <  c  id 

Lot  u»  partition  C  into  kd  interval*  of  oqual  l*ngth> 

Ej  »  (  8  I  i*k~d  <  6  <  (i+l)*k"d  )>  i-0,1,  ...  ,  kd-l 

If  th*  distune*  b*tw**n  any  two  con**cutiv*  point*  alony  tha 
h*uri*tic  tour  i»  >  f<k"d><  th*n,  by  <P2>,  th***  point* 
cannot  li*  within  th*  »am*  Ej>  »o 

#(X(A  I  v  >  f<k“d>>  <  #<i  I  <*<£*>  n  S>  i*  non*npty>  <6.2> 

But  th*  point*  in  any  E*  li*  within  f<k”2>  of  *ach 


oth*r,  *o 


•a 
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<*KE*[>  n  S)  it  nonempty  -fr  tCE^)  c  S<f<k“d>>.  <6. 3) 

W  «H 

Finally >  aach  E^  haa  langth  k  >  and  i*  maaaura 
proaarving;  *o 

VC«KE![>3  «  k"d  <8.  4> 

Containing  <8.0  through  <8.4),  and  aatting  t  *  f<k_d) 

*  ft/' k,  wa  obtain 

#<X«A  I  x  >  ft^k>  <  kd  VCS<ft/k)3  «  c *  ft*”1  k  L*  +  c  ftd 

If  £  it  i  nonnag at iva  nunbar  <not  nacavtarily  intagar),  than 

#<XtA  l  x  >  ft/k>  <  c'  ft61"1  f k |  L*  «•  c  ftd 

<  c'  ftd'1  <k+l>  L*  +  c  ftd 

Lotting  t  ■  ft/k,  this  baeonaa 

#<X<A  I  x>t>  <  ftd  <c'  L*Vt  +  c'  L*'(l  +  c> 

i 
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There  are  N  points  in  S>  so 

#<X«A  I  X>t>  <  N. 

The  distinct  between  any  two  points  in  S  is  bounded  above  by 
f(l/2>  (since  1/2  maximizes  f(O)  and  by  L*/2  (because  the 
triangle  Inequality  holds  and  an  optimal  tour  of  length  L* 
joins  the  points  twice).  Thus 

#<x«A  l  x>t>  -  0,  t  >  min(f>  2'6,  L*/2>. 


These  combine  to  form 


‘  JC 


min(rt  2~d,  L*'2> 


So 


L/L%  -  JfL  2  min(Nxl_*> 


■  c 


,  Ad  r 1^2ftdc' 


z 


min(N>  Ad  <c'L*>t  +  c'L*/ft  +  c>)  dt 

c'Ad't>  dt  +  c'  <.Ct'2>  6  +  cfid/2 

min(N/l'/T>  dx  +  constant  terms 

<t  -  t^c'ftdL* 


> 


■  c'  ftd  ln(N)  +  constant  terms.  D 


(6.5) 
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?,  Proof  of  Thtortn  3.1.  We  consider  first  a 
spacefilling  curve  closely  related  to  <♦»*  and  defined  by  a 
similar  recursive  structure  in  which  the  domain  is  broken  into 
two*  rather  than  four*  identical  parts.  This  related 
spacefilling  function  will  be  denoted  by  its  domain  is 
the  set 


D  *  <  w  I  0  <  «  <  2  >;  <7.1 > 

and  its  range  is  the  triangle 

T  ■  <  (x*y)  I  x+y<2j  x*y>0  ).  (7.2) 

fls  «  ranges  from  0  to  2*  the  related  spacefilling  curve 
<k«)>  will  continuously  cover  every  point  in  T*  starting  at 
the  lower  right  vertex  (2*0)  and  ending  at  the  upper  left  vertex 
(0*2). 

We  define  *  recursively*  by  splitting  T  into  two 
smaller  triangles*  each  filled  with  a  spacefilling  curve 
(Figure  3).  The  smaller  triangles  are  given  by 


Tq  ■  <  <x*y)«T  I  x>y  ) 
Tj  *  <  <x*y)iT  I  x<y  > 


(7.3) 
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Me  Also  divide  the  domain  D  into  two  helve* 


Dg  ■  (  ««D  l  «< l  > 

D,  *  <  uiD  I  «>1  > 
1  — 


C?.4> 


»o  that  ♦  MAp*  each  onto  T..  Finally,  we  construct 
Affine  Mapping*  Mi  froM  T  to  and  from  Di  to  D, 
each  conpletely  determined  by  the  following  information 


R0  •  <2,0>  <2,0>,  C0,0>  +  Cl,l>,  <0,2)  +  C0,0> 

Rj  ••  <2,0>  ♦  C0,0>,  <0;0>  +  Cl,l>,  CO, 2)  *  CO , 2)  C7.3> 

b0  ■  0  ■>  0,  1  ■>  2  b*  •  1  ♦  0,  2  ♦  2 

The  recur*lve  *trueture  of  ♦  may  now  be  written  a* 


i*0, 1 


C7.6> 
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Lemma  7. 1.  Eq.  0.6?  has  a  unique  solution  t  if  and  only 
if  <7. 6)  has  a  unique  solution  4>.  Furthermore,  if  both 
solutions  ax 1st,  than 

♦  <8>  *  <1,1?  -  *<6+.5>,  8«C  <7.7? 

Proof .  Lat  8  ■  <J+h?r8,  where  j  is  an  integer.,  and 
0<h<l.  Each  case  J*0 ,  1 , . . . ,  7  Must  ba  considered  separately. 

Me  sketch  the  case  j» 0  only.  Suppose  first  that  ■+■  exists 
and  <7. 7?  holds.  The  LHS  of  <3.6?  is  'KB  J-'Kh/’S? .  By  <7.7>, 
this  bacOMes  <  1  >  1  >  -  ♦<.5+h/'8).  But  the  argument  of  * 
is  lass  than  one,  so  <7.  £3  permits  this  expression  to  be 
converted  to  <1,1?  -  fig £■$>< l+h/4? 1 .  Continuing  in  this 
Mannar.,  the  LHS  of  <3.6?  is  seen  to  ba  <1.1?  - 
fl0Cfl1rfl0t«0C'V<2h?3]3].  The  RHS  of  <3. 8) 
becomes  RgCCi.  i?-fl1cfl0C'f<2h>333.  These 

expressions  are  equivalent.  Similarly,  if  we  assume  that  * 
exists,  we  can  construct  ♦  .  □ 

Eq.  <7.6?  may  also  be  written  as  a  fixed-point  identity 

•a  ■  F  4>  <7.8? 


where  F  is  a  transformation  on  the  space  of  T -valued  functions 
on  D.  For  any  function  ♦  •  D  ■»  T  ,  the  image  of  ♦  under  F 
is  given  by 
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CF43<»)  -  fiiC4<biC»3)3>  i-0,1  <?.9> 

The  operator  F  provide*  *  basis  for  demonstrating  the  existence 
of  *. 


Lenn*  7.2.  The  trfl.n5forn4t.ion5  and  bj  fl.re  distfl.nce 
preserving,  in  the  following  sense; 

a)  II  fliC<x,y>3  -  FliC<xSyi  >3II2  ■  2“1/2  IKx,y)-<x',y '  )IU, 

for  any  <x,y>, <x' >y' >«T,  and  i=0  or  1. 

b?  Ib^Cal  -  b.C«'.7l  »  2  lo-o'l, 

for  any  and  1*0  or  1. 

Corollary  7.2R,  F  is  a  contraction  operator  under  the  *up 
<Euclidean>  nom,  i.e., 

5UpWfD<IICF4Ko)-CF4/3Cw)|l2>  <  2~1/2supM(D<  !!♦<«>-♦' Co) flgl . 

Coroll aru  7.2B.  For  any  "initial  guess"  ♦,D+T,  the 
sequence  of  function#  n*0,l,  ...  ,  converge#  unifornly 

to  a  (uniquei  solution  of  <7.8>  tequivalantly >  of  <7.6)3. 
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Corollaru  7.2C.  There  exists  *  unique  continuous  function  f 
satisfying  (3.6). 

Proof 1  Lot  ♦  be  th®  unique  continuous  solution  of  <7. 8) 
and  apply  L®nna  7.1.  □ 


fl  siMple  initial  guess  for  t  is 

♦  <m>  •»  «)>  o«D  <7.io) 

Lot  »  Fn4.  By  Corollary  7.ZB,  4.  Figure 

4  show*  th®  first  f®w  t®rn*  of  this  *®qu®nc®.  In  vi®w  of  <7. 7). 
th®  r®striction  of  to  S  (shown  within  dotted  lin®*  in 
Figure  4)  is  a  convergent  sequence  of  approxi nation*  to 
rotated  by  180°. 


to  the  spacefilling  curve  \Jj .  The  approximations  to  ij  are  shown 
within  the  dotted  lines. 
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Me  now  turn  our  attention  to  establishing  <P1>. 

Lemma  7,3,  If  X  and  Y  are  tooth  integer  Multiple*  of 
2-k ,  then 

<a>  *<THETfl<X,Y>)  -  <X, Y> . 

<b>  TH£Tfl<X,Y)  will  toe  an  integer  multiple  of  2~2<k*i:,t 
<c>  THETfl<X,Y>  will  call  it*elf  at  most  k  tine*. 

Proof 1  Pin  induction  on  k  follow*  easily  from  the  structure 
of  R.  and  0 

Corollary  7. 3f)  The  function  ♦  *ati*fie*  property  <P1>. 

Corollary  7.3B.  Y  is  a  spacefilling  curve  over  S. 

Proof  =  Let  §  denote  the  point*  in  S  who*e  coordinate* 

have  finite  binary  expan*lon.  THETfl  will  compute,  for 

<x,y)«$,  a  8  *uch  that  Y<8>  ■  <x,y>.  Consequently,  § 

c  <KC)  c  S.  C  i*  a  compact  set,  so  the  continuity 

of  Y  <a**ured  toy  Corollary  7.20  implies  that  *KC)  is  compact. 

That  is,  closure<§>  c  tCC>.  But  closured)  =  S, 

so  t<0  -  S.  Since  Y  is  continuous  and  surjective,  it  is 

spacefilling.  0 
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Next/  we  establish  <P2>. 

lemma  7.4.  >*2  <  2  i  to-»'  I  *  «#/o'€D. 

Proof.  Since  V"**,  it  suffices  to  show  that 

>«2  <  2  J  I «-« *  1 ,  6>,6)*tO,  n“0  si-,.,  <7. 1 1  > 

This  we  do  by  induction.  By  <7.10)> 

!l*0<e>-*0<e'>l»o  <  JT  •  iw-«'l/  w/»'«D. 

tm  — 

But  <7.1)  ensures  that  lo-w'  I  <  2/  W/«'«D/  and  so 

<7. 11>  holds  when  n=0.  To  perform  the  induction/  we  consider 

two  case*1 

<a>  «/«'<l  or  «/«'>!.  The  induction  follows 

am 

trivially  from  Lemma  7.2. 

<b>  »<i<«'  Of  «'<1<«.  Since  any  two 
points  in  T  determine  an  angle  at  <0,0  >  of/  at  most/  90°/  the 
"cosine  law"  extension  of  the  Pythagorean  Theorem  implies 
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>II2]Z 

£  <o>-<0  >0  >  llg]  2  ♦  [fKO.OJ-f^^w'JHgj2 

But  4r,+1<t>-<0,0>  [«•«  <?.  10>,  <7. 11>  and  <7.5>,  or 
axaMina  Flyura  33  amd  «o  this  nay  ba  written  a* 

(ll4n‘flC«>-41f>+1<ttOII2j2 

<  ^M4T,4l<»)-4T'+1<l>ll2]2  ♦  [»41n4l<l>-4in+1<aOII 

FroM  <a)>  above*  it  follow*  that 

2 

(ll4T,4^<»>-4tvf*<tt'  >llg]  <  4*  ltt-1 1  +4-I1-WM  -4-lw 


Thi*  establishes  the  induction  for  case  <b>.  0 


Corollary  7,4ft  The  function  *  satisfies  <P2>  with  11*11  * 
Euclidean  distance  in  <2.1>  and  ft  «  2  in  <2.2>. 

How  we  establish  <P3>. 

Lemma  7.S,  For  any  integers  m>t  and  0<i<2M*  the 
set  <  *<•>  I  <  t  <  <i+l>*2“M  >  is  a  right  triangle* 

of  area  2‘’‘^^  having  45°  angles  at  ♦<i*2~M>  and 
*<<i+l>*2"ri>*  and  a  90°  angle  at  'K<i+.5>-2""rn). 

Proof 1  The  assertion  is  true  of  t  for  m«0  by  definition 
of  T.  fln  induction  on  m  is  based  on  the  recursive  structure 
shown  in  Figure  3«  as  m«-m+l>  each  triangle  is  bisected  at  its 
right  angle  to  form  two  new  triangles  of  equal  area.  In  view  of 
Lemma  1*  the  property  extends  to  *  as  well.  □ 

Corollary  7.3fl.  t  is  Lebesgue  measure  preserving. 
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Finally/  we  examine  the  sequence  of  meant/ 

Lemma  7.6,  For  any  N/  N' / 

1  i  ^  I  lnCN-'N' >  I . 

Proof  >  If  a  point  selected  at  r*ndon  it  removed  from  a 
heuristic  tour  of  N  independent  random  points/  the  remaining  N-l 
points  will  be  independent  and  sequenced  so  as  to  form  a 
heuristic  tour.  By  the  triangle  inequality/  the  heuristic  tour 
of  N-l  points  is  no  longer  than  the  heuristic  tour  of  N  points/ 
but  it  is  less  by/  at  most/  the  lengths  of  the  two  links 
adjoining  the  removed  point.  8y  symmetry/  the  expected  length  of 
a  link  selected  at  random  it  N'1  times  the  tour  length. 

That  is* 

0  <  N1'2  pN  -  <N-1)1/2  ln-1  <  <2rN>  N1'2  pN 

Dividing  by  N*/2/  and  using  1-N"1  <  <  1/ 

as  well  as  <3. 5>/  this  becomes 

I^N  ~  ^H-l1  - 
For  N'<N/  then/ 

'“H  “  **H'  1  i  Zi-N'+l  i0/i  i  2“'*  dt  D 
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Lemma  ?.  ?.  There  i*  «  sequence  such  that 
/*4N  *  for  Al1  N.  and  p^-p^  +  0. 

Proof »  Consider  a  Modification  of  the  random  N  point 
problem.  where  the  number  of  point* .  denoted  by  n.  1*  it** If  a 
random  variable  having  Poi««on  di*trtbution  with  mean  N.  The 
8^'s  now  constitute  a  Poi»*on  process  on  C.  Let  p' <N>  * 

1  y  Q 

H  •  the  expected  1 anyth  of  the  heuristic  tour  of 

these  n  point*.  Also  lot  p*(N)  *  N”1^2  ■  the  exp  act  ad 

length  of  a  heuristic  tour  of  the  n  point*  along  with  the  four 

vertices  of  the  square  (that  is.-  n+4  points  in  all).  and  let 

n"<N)  *  N”1^2  •  the  expected  length  of  all  link*  along 

the  heuri*tic  tour  of  the  n  randoM  points  and  four  vertices 

excepting  tho*e  link*  that  adjoin  a  vertex  of  the  square.  Now 

0  <  p+(N)  -  p"”<N)  <  since  <i>  the  tour* 

described  by  p*(N>  and  p“(N)  differ  by  (at  Most)  the 

length  of  eight  link*  (adjoining  the  vertices  of  the  square). 

each  of  which  range*  over  1'N  of  C  on  average,  and  <ii>  f<0  i* 

concave,  so  that  ECf<4)3  <  f(ECA3>.  The  recursive 

structure  of  4  <*hown  in  Figure  l>  assure*  that  p”(N>  < 

p“<4N>  and  p+<N)  >  p+(4N).  so  p^N^e1*)  and 

A  I* 

p  (N*4  )  converge  nonotonically  (as  k+»)  to  a  limit 
(that  nay  depend  on  N).  Let  pN  denote  this  limit.  Clearly 
p4N  ■  pN.  *lso  p""(N)  <  <  p+(N>.  find, 

by  their  definition  (above).  p“<N)  <  p'(N)  < 
p+(N).  80 
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Ipt' <H>  -  <  H+<N>V<N>  <  ®G/N« 

The  ralAtionship  between  <N>  end  HN  i* 

^<M>  -  EC<ti/N)1/£  *  Et<r»/N>1'"2  V-^ 

.  £CCnfN>1'8  <1*„-HN>3  *  «<«/t4>1/2  *  '=  t*M  *  ^ 


Thi#>  Along  with  LemmA  7.6  And  <2.5>  implies 

m- <H>-t^,t  <  20  Et<«>/N)1/2  Ilncn/N>i:  +  ft  ECI<n/N>‘'2  -  1H. 
Wh.„  N  1.  U n>.  "'N  U  normally  dl.tributad  with  "•*"*>  «* 

varlaneaalfN,  »«i  *°  D 


mnnllaru  7.7ft.  Thar.  i»  a  conti nuou.  function  .<■>  *“ct1 


thAt 


^log4<N)  mod 


+  0. 


>  log4<N>  mod  is 

Proof !  Let  pl4  *  *  V 

Continuity  of  p<->  follows  from  LemmA  7.6.  □ 


The  various 


Assertions  of  Theorem  3. 1  have  now  been  proved. 
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6.  Conclusion*.  The  spacefilling  heuristic  should  prove 
useful  in  large  applications  because  of  its  speeds  only 
CKN  log  N)  operations*  an  order  of  Magnitude  faster  than  any 
other  T3F  heuristic  coMnonly  considered.  It  achieves  this  by  the 
surprising  tactic  of  ignoring  the  0<N2*>  interpoint  distances. 

Our  Method  is  also  robust*  in  the  sense  that  it  constructs 
a  single  tour  which  is  good  with  respect  to  a  variety  of  Metrics 
and  point  distributions.  Thus*  problem  parameters  required  by 
other  Methods  need  not  be  measured*  and  need  not  remain  constant 
over  time.  Furthermore*  since  it  is  based  on  sorting*  our 
heuristic  has  the  advantage  that  insertions  and  deletions  can  be 
made  in  only  OClog  N>  operations. 

The  siMplicity  of  our  heuristic  should  make  it  attractive 
for  low  technology  applications.  It  requires  no  real 
Multiplications  or  square  roots*  and  so  should  execute  quickly 
on  microprocessor-based  systems.  In  fact*  a  coMputer  isn't  even 
required.  For  example*  a  table  o**  precomputed  0  values  can  be 
used  to  convert  <x*y>  coordinates  (read  from  a  map)  to  numbers 
that  are  then  sorted  manually.  With  our  students  Lee  Collins  and 
Gill  Warden*  we  di vised  just  such  a  routing  system  for  an 
Atlanta  charitable  organization.  It  delivers  meals  to  over  200 
locations  which  change  daily.  Our  heuristic  is  at  the  heart  of 
this  system*  and  is  implemented  on  two  card  files.  The  system  is 
easily  maintained*  cost  virtually  nothing*  yet  reduced  total 
mileage  by  \3>,  and  enabled  a  reduction  from  five  to  four 
delivery  vehicles.  For  details*  see  143. 
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APPENDIX 

The  following  BASIC  code  will  compute  a  tour*  of  N  point*  in 
a  square  of  side  S.  Variables  whose  names  start  with  I-N  will 
contain  integer  values.  This  program  may  require  modif ications 
to  run  under  some  versions  of  BASIC  <e.g.  eliminate  variable 
dimensions  in  30  and  410 >.  Our  sorting  algorithm  is  HEAPSORT 
<see  Kmthttl]). 


The  structure  of  the  program  is  as  follows1 


Statements  10-90 
Statements  100-150 
Statements  200-330 
Statements  400-820 


input  and  output  of  variables 
calculation  of  for  <x^y^> 
sort  Oj's 

optional  improvement  procedure 
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10  INPUT  "NUMBER  OF  POINTS";  N  ■  K»10  '  KP  -  2t<K-l> 
20  DIM  AfCN),XCN>,YCN),THCN),  IQOO,NRCN> 

30  INPUT  "SIDE  OF  SQUARE";  S  >  U-S*.501/KP 
40  PRINT  "NAME,  X,  Y«"  1  FOR  I-I  TO  N 

.  INPUT  Af<I),X<n,Y<n  •  NEXT  I 
50  FOR  1-1  TO  N  ■  GOSUB  100  •  NRC I  )«I  *  NEXT  I 
S®  GOSUB  200  *  GOSUB  400 

?0  PRINT  -  PRINT  " RANK "  > " NAME "  ,  "  X  " ,  "  Y  " , " THETA " 

80  FOR  1-1  TO  N  :  J»NR< I > 

«  PRINT  I  j  Af  <  J  > ,  XC  J ) >  YC  J  >  >  THC  J  )  *  NEXT  I 

90  STOP 


100  REM  *#*  SUBROUTINE  TO  COMPUTE  THC I >  *## 

110  KX-INTCXCI  VU>  *  KY-INTCYC  I  VU> 

120  FOR  J-l  TO  K  ■  JX-INTCKX/KP?  <  JY-INTC KY/KP ) 

«  KX-2*C KX-KP* JX )  -  KY«2*< KY-KP* JY > 

122  *  IQC J  >JY+0*  JX-2#  JX*  JY  >  NEXT  J 

130  T-IQCKV4 

140  FOR  J*K-1  TO  1  STEP  -1  •  T*T+C 6-IQC J )  V4 

*  T-T-INTCT)  >  T*C  3. 5+T+IQC  J )  V4  ■  NEXT  J 
150  THC  IW-INTCT)  «  RETURN 


200  REM  ***  SUBROUTINE  TO  SORT  ARRAY  NR  SO  THAT  *** 

210  REM  ***  THC  NRC  1  >  X  THC  NRC  2  >  >< .  .  .  <  THC  NRC  N  >  >  *** 

220  IL*INTCN/2)+l  >  IR-N 

230  IF  IL>1  THEN  IL-IL-1  >  NA«NRCIL>  !  GOTO  260 
240  NA-NRCIR>  *  NRCIR>-NRC1>  '  IR-IR-1 
250  IF  IR«1  THEN  NRC  15-NA  •  RETURN 
260  TA-THCNA)  ■  J-IL 

270  I-J  -  J-2*J  '  IF  J«IR  THEN  GOTO  300 
280  IF  J>IR  THEN  GOTO  320 
2 90  IF  THC  NRC  J  >  X THC  NRC  J+ 1 >  >  THEN  J-J+l 
300  IF  TA>»THC NRC J))  THEN  GOTO  320 
310  NRC I )»NRC  J )  <  GOTO  270 
320  NRCI)-NA  «  GOTO  230 
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400  REM  ***  OPTIONAL  IMPROVEMENT  SECTION  *## 

410  DIM  LF<N>,LB<N?,D2<3?,KV<3?,KW<3>  •  KR«KP*2-1 
420  NR(  0  >«NR<  N  )  •  FOR  1  =  1  TO  N  >  LF< NR< 1-1 > >=NR< I > 

:  LSCNRC I >?*NR< 1-1 >  »  NEXT  I  =  REM  ***  SET  UP  LINKS  *** 
430  FOR  11  =  1  TO  N  >  JJ=NR<  II)  =  MX*X<JJVU  ■  MY=Y<JJ?/U 
440  T=TH(  JJ ?  :  REM  *#*  GET  TRIANGLE  CONTAINING  POINT  II  *** 

450  TB=TH< LB< J J > ?  «  IF  TB>T  THEN  TB=0 

460  TF-THC LF< JJ > >  =  IF  TF<T  THEN  TF=1 

470  TL--.5  «  TG-1.5  =  KV< 1 ?=-KR  «  KWC  1  >=KR 

'  KV(  2  )»KR  »  KM<  2  >=KR  =  KV<3?=KR  =  KW<3?»-KR 
480  TM=< TL+TG )/2  «  IF  TB<TL  AND  TF>TG  THEN  GOTO  560 

490  IF  T>TM  THEN  GOTO  510 
500  TG=TM  «  IV=3  •  GOTO  520 

510  TL-TM  :  IV«1 

520  KK=<  KV<  1  ?+KV<  3  )  >/2  <  KV<  IV  ?=KV<  2  >  •  KVC  2  ?=KK 

530  KK»<  KW<  1  >+KW<  3  ?  '>/2  ■  KW<  IV  ?=KW<  2  >  •  KW<2?=KK 

540  GOTO  480 

550  REM  ***  TRIANGLE  VERTICES  ARE  CU*KVC .  >,U*KW< .  ?>  *** 

560  FOR  IV-1  TO  3  «  D2<  IV?*<MX~KV<  IV>)f2+CMY-KW(IV>>t2 

!  NEXT  IV  >  KV<  0  ?*KV<  3 )  •  KW<0?=KW<3?  «  D2<0?=D2<3> 

570  DS=1E9  «  IS=0  :  REM  **#  TEST  EACH  PROJ'N  *** 

580  FOR  IV«1  TO  3 

590  RD=<  <  KV<  IV  )-KV<  IV-1  ))f2+<KW<  IV  >-KW<  IV-1  >>f2  ) 

600  R1»C1-<D2<  IV>-D2C  IV-1  ?  VRDV2  '  R2-1-R1 

610  KX*R  1  *KV<  I V  ?+R2*KV<  IV-1?  «  XC  0  >=U*<  KX+SGNC  KX-MX  ?  ? 

620  KY=R1*KW< IV?+R2*KW(  IV-1 ?  >  Y< 0 >=U*< KY+SGNC KY-MY > > 

630  1=0  •  GOSUB  100  *  REM  ***  FIND  THETA  FOR  PROJ'N  *** 

640  T=TH<0?  ■  IF  T>TB  AND  T<TF  THEN  730 

650  NL*1  >  NG=N  =  TL«0  =  TG=1 

660  NM=INT<  <  NL+NG  )/2  ?  «  TM«TH< NR< NM > >  •  IF  NG-NL<2  THEN  690 
670  IF  T>TM  THEN  TL-TM  >  NL=NM  «  GOTO  660 

680  TG«TM  :  NG=NM  *  GOTO  660 

690  JL»NR<NL?  !  JG-NRCNG?  ■■  REM  ***  NEW  POS'N  TO  TEST  *** 

700  IF  LF<  JL  ?<  > JG  OR  JL-JJ  OR  JG=JJ  THEN  GOTO  730 

710  O-SQRC  <  X<  JL  >-X<  JJ  )  )t 2+<  Y<  JL  >~Y<  JJ  >  >t2  > 

+SQR<  <  X(  JG  ?-X<  JJ  )  ?  t2+<  Y<  JG  >-Y<  J  J  >  >*2  ? 

-SQRC  <  XC  JG  >-X<  JL  ?  >t2-K  Y<  JG  ?-Y<  JL  ?  >t2  ? 

720  IF  D<DS  THEN  DS=D  >  IS=IV  =  TS»T  <  JS-JL 

730  NEXT  IV 

740  IF  IS»0  THEN  GOTO  800 

750  JB=LB<  JJ  >  ■  JF=LF<  JJ  >  >  REM  **#  BETTER  THAN  OLD  POS'N?  *** 

760  D=SQR<  < X<  JF  ?-X( JJ ? }t2+<  Y< JF  >-Y<  JJ  ?  ?t2  ? 

+SQR<  <  X<  JB  >-X<  JJ  )  ?t2+<  Y<  JB  >-Y<  JJ  >  >72  > 

-SQR<  <  X<  JF  >-Xf  JB  ?  >t2-K  Y<  JF  >-Y<  JB  )  ?t2  ? 

770  IF  D<=DS  THEN  GOTO  800 

780  NR<  1 1  >*LF<  JJ  ?  :  TH<  J  J  )=TS  =  JH=LF<JS>  >  LB<JJ?*JS 

790  LF<JJ?=JH  =  LF<  JS  )=JJ  *  LB<JH?=JJ  «  LF<JB?«JF  =  LB<  JF  )*JB 

800  NEXT  I I 

810  J=NR< 1 )  <  FOR  1=2  TO  N  :  J=LF<J>  =  NR< I ?«J  :  NEXT  I 

820  RETURN 


